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Abstract For a sample of metal-poor stars (—3.3 < [Fe/H] < —2.2) that have high- 
resolution spectroscopic abundance determinations, we have measured equivalent widths 
(EW) of the Ca II K, Mg I b and near-infrared (NIR) Ca II triplet lines using low-resolution 
spectra of the Sloan Digital Sky Survey (SDSS), calculated effective temperatures from (g — 
z)o color, deduced stellar surface gravities by fitting stellar isochrones, and determined metal- 
licities based on the aforementioned quantities. Metallicities thus derived from the Ca II K 
line are in much better agreement with the results determined from high-resolution spectra 
than the values given in the SDSS Data Release 7 (DR7). The metallicities derived from the 
Mg I b lines have a large dispersion owing to the large measurement errors, whereas those 
deduced from the Ca II triplet lines are too high due to both non-local thermodynamic al equi- 
librium (NLTE) effects and measurement errors. Abundances after corrected for the NLTE 
effect for the Mg I b lines and Ca II triplet lines are also presented. Following this method, 
we have identified six candidates of ultra-metal-poor stars with [Fe/H] ~ —4.0 from a sam- 
ple of 166 metal-poor star candidates. One of them, SDSS J102915+172927, was recently 
confirmed to be an ultra-metal-poor ([Fe/H] < —4.0) star with the lowest metallicity ever 
measured. Follow-up high-resolution spectroscopy for the other five ultra-metal-poor stars in 
our sample will therefore be of great interest. 
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1 INTRODUCTION 

Elemental abundance patterns in the cosmos result from the nuclear processes that happened in the Big 



Bang, stars, and the interstellar medium (Pagel 1997 1. Measurements and analyses of stellar element abun- 



dances can help to address the formation and evolution of stars and galaxies. Searching for extremely 
metal -poor (EMP, [Fe/H] < — 3; Beers & Christlieb 2005) stars and follow-up detailed analyses of their 
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elemental abundances can give us a deep insight of the formation of the first generation of stars and 
the early chemical evolution of galaxies (Frebel 120101) . To obtain accurate elemental abundance of EMP 
stars, one needs time-consuming high-resolution spectroscopic observations. As such, candidates of EMP 
stars are first selected from low to intermediate-resolution spectra for follow-up high-resolution spectro- 
scopic observations (Beers & Christlieb 12005b . The most iron-poor stars discovered so far are all candi- 
dates selected by low-resolution spectroscopic surveys, such as the HK survey JBeers et al. 19921 1 and the 
Hamburg/ESO survey jWisotzki et al. 19961 , and their ultra-low metallicities are later confirmed by follow- 
up high-resolution spectroscopic observations, e.g., HE 0107-5240 ([Fe/H] = — 5.2; lChristlieb et al. 2002t . 
HE 1327-2326 ([Fe/H] = -5.4: IFrebel et al. 2005llAoki et al. 20061 . and HE 0557-4840 ([Fe/H] = -4.75; 
INorris et al. 2007b . Thus, reliable metallicity determinations for metal-poor stars from low to intermediate- 
resolution spectra is the key to find and identify more EMP stars. 

Using spectra of a resolving power R ~ 7500 from the Radial Velocity Experiment (RAVE) survey, 
Fulbright et al. (12010b develop an efficient technique of identifying candidates of ultra-metal-poor stars 
([Fe/H] < —4) from a pool of metal-poor ([Fe/H] < —2) stars, without the need of further observations of 
higher spectral resolutions. With this method, three stars with [Fe/H] < —4 were identified in a sample of 
approximately 200,000 stars. 

The large number of low-resolution spectra provided by the SDSS and to be provided by the forthcoming 
LAMOST (The Large Sky Area Multi-Object Fiber Spectroscopic Telescope, also named The Guoshoujing 
Telescope; |Wang et al. 1996] ISu et al. 19981 |Xing et al. 19981 IZhao 20001 ICui et al. 20041 IZhu et al. 20061 



IZhao et al. 20121 c.f. http://www.lamost.org/website/en/) Galactic surveys (Deng et al. 2012 Liu et al. 2012, 
in preparation) motivate the need for better metallicity determinations from those spectra, in particular at 
low metallicity, such that reliable stellar candidates of EMP stars can be identified for follow-up high spec- 
tral resolution, high signal-to-noise ratio (S/N) spectroscopy. Preliminary observations of LAMOST ob- 
tained during the early commissioning phase have shown its ability of searching for metal-poor stars in 
the Milky Way dLi et al. 20101 IWu et al. 20101 iRen et al. 2012b . As one of the hitherto most successful and 
productive surveys, the SDSS has accumulated hundreds of thousands low resolution stellar spectra that 
can be utilized to search for candidates of EMP stars. The stellar atmosphere parameters released in DR7 
of the SDSS ( |Abazajian et al. 2009b, are provided by the SEGUE Stellar Parameter Pipeline (SSPP). For 
the DR7 version of the SSPP, the metallicity determinations at the extremum of [Fe/H] < —3.0 were not 
well constrained due to lack of corresponding calibrators (Smolinski et al. 1201 lb . There is evidence that 
the SSPP-derived metallicities for stars below [Fe/H] = —2.7 can be 0.3 dex higher than the actual values 



determined by high-resolution spectroscopy (Yanny et al. 2009 1. The performance of the SSPP has been 



significantly improved in DR8 of the SDSS (lAihara et al. 201 ft . yet the most metal deficient star found by 
the SSPP in this huge pool of more than half a million stars has an unexpected, rather high metallicity, 
[Fe/H] = —3.9, indicating that even the current, most up-to-date version of the SSPP might have missed the 
most metal-poor stars. Motivated by the success of Fulbright et al. (20101 in finding EMPs with the RAVE 
spectra, we have applied their technique to the huge data set of the SDSS spectra of even lower spectral 
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resolution than those of RAVE. We show that the method is not only applicable to the SDSS spectra, it is 
even capable of finding candidates of [Fe/H] < —4. 

The paper is organized as follows. In Section 2, we describe our method and apply it to a sample of 13 
metal-poor stars that have both the SDSS and high-resolution spectroscopic observations. In Section 3, we 
extend the analysis to a large sample of low-metallicity star candidates in the SDSS data archive and search 
for candidates of EMP stars. Our summary follows in Section 4. 

2 METHOD 

2.1 Selection of Sample Metal-Poor Stars with High-Resolution Spectra 

The abundance data derived from high-resolution spectroscopy of known metal-poor stars can be found in 
the table that contains elemental abundances of ~ 1000 metal-poor stars collected from the literature by 
Frebel (1201 Ok and the Stellar Abundances for Galactic Archeology (SAGA) Database dSuda et al. 20081) 
that contains 1212 metal-poor halo stars with [Fe/H] < —2.5. Those valuable data provide an effective 
way to examine our method and test the results. We select metal-poor stars with [Fe/H] < —2.0 from the 
compilation of Frebel and the SAGA database, and identify them in the SDSS archive. In total, we find 13 
metal-poor stars with the SDSS spectra. Table [D lists the stellar parameters deduced from high-resolution 
spectroscopic analyses and those given in SDSS DR8. Note that there exist systematic uncertainties in the 
high-resolution abundance data collected from different literature studies (the original references are also 
listed in TableHJ) (Frebel l2010l ). But for the relatively small-scale systematic uncertainties (around 0.3 dex, 
Frebel 1201 01 ) will not affect our conclusions, we can use this high-resolution reference sample to provide 
robust estimates for [Fe/H]. 

A comparison of the abundances derived from high-resolution spectra and those determined with the SSPP 
used in DR8, it is clear that the SSPP has overestimated the metallicities of those metal-poor stars, by ~ 
0.1 dex on average (see also Fig. |3(d)| in Section l23l ). In addition, the SSPP also seems to have overestimated 
the effective temperatures (by approximately 100 K) than those adopted in high-resolution spectroscopic 
analyses (c.f. Fig. |l(b)| in Section 12.31 ). The latter is probably an important factor leading to the higher 
metallicities of the SSPP. 

2.2 The Adopted Metal Line Features 

Given the limited S/N and spectral resolution of the SDSS spectra, weak Fe I lines can only be used for 
metallicity analysis in a limited stellar parameter range dRen et al. 20 12). The Ca II K line at 3933.7 A, the 
Mglb lines at 5167.3, 5172.7 and 5183.6 A, and the Call triplet lines at 8498.0, 8542.1 and 8662.1 A, are 
on other hand stronger and easier to measure for metal-poor stars. The Ca II K line is probably the best indi- 
cator for the overall metallicity of a metal-poor star when only low S/N and resolution spectra are available 
(Frebel 1201 01 ). Generally, there exists a degeneracy amongst the stellar parameters including effective tem- 
perature, surface gravity and metallicity. A lower-metallicity, cooler giant and a higher-metallicity hotter 
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Table 1: Stellar Parameters of the Sample of 13 Stars. 
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turnoff star can have equal strengths of those lines (Fulbright et al. 2010). For the Ca II triplet lines, a 100 K 
difference in temperature leads to about 0.1 dex difference in metallicity (Fulbright et al. ,2010), i.e., an 
overestimated effective temperature will result in an overestimated metallicity. 

For metal-poor star of [Fe/H] ~ —1.5, one can assume that the cv-element (Ca, Mg, etc.) abundances are 
enhanced by 0.4 dex with respect to the Solar cv/Fe ratio (Fulbright et al. 120101 1. Although there are some 
outliers that have a strong enhancement of a element (Frebel 120101 fig. 2), this assumption is generally 
found to be valid in most cases and is also assumed in the current work when we convert the deduced 
calcium or magnesium abundance to that of iron. 

2.3 Determination of Stellar Atmosphere Parameters 

2.3.1 Effective Temperature 

Being the most important factor influencing estimates of metallicities, erroneously assigned effective tem- 
peratures can result in biased elemental abundances. For our sample, the DR8 version of the SSPP pro- 
duces a higher estimate of effective temperature, ~ 100 K higher than those adopted in high-resolution 
spectroscopic analyses. Fulbright et al. (2010) utilize the NIR photometry from the 2MASS survey 
(ICutri et al. 20031 1 and optical UBV photometry to calculate photometric effective temperature r p hot- Here 
we make use of the high quality SDSS ugriz photometry and use the long-baseline (g — z )o color as effec- 
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tive temperature indicator. Two calibrations are available to us. 

Ludwig et al. (2008) obtain a theoretical relation between the (g — z)o color and the effective temperature 
based on model atmospheres, 

r ofr = 7126.7 - 2844.2 • (g - z)o + 666.80 • (g - z)l - 11.724 • (g - z)l (1) 

Smolinski et al. ( 1201 II ). based on likely cluster members and calibration, derive two color-temperature 
relations applicable to different metallicity regions. For the range of elemental abundance of interest here, 
we adopt their relation for [Fe/H] < —1.5, 

T cS = 6993 - 2573 • (g - z)o + 530.9 • (g ~ z) 2 . (2) 

The effective temperatures derived from the above two calibrations are referred to as T\ and T?. T\ and 
T2 are respectively 31 and 74 K lower than the values adopted in high-resolution spectroscopic analyses 
(mostly are photometric temperatures) on average. 

2.3.2 Surface Gravity 



We adopt the Yonsei-Yale isochrones (Demarque et al. 2004 1 for the 13 metal-poor stars. Considering that 



most metal-poor stars are old stars, we adopt an age of 12 Gyr for all of our stars. As the metallicity can- 
not be easily derived by measuring the equivalent width of faint Fe I lines in the SDSS spectra, we use 
the metallicity estimates in SDSS DR8 as initial values. The values of surface gravity obtained from the 
isochrone assuming the effective temperature calculated above are then set as the initial values of stellar 
surface gravity (see Fig. [2j. Guided by results adopted in the high-resolution spectroscopic analyses, for 
cooler stars with log T c s < 3.79, we use the upper parts of the isochrones (above the turn-off points), while 
for hotter stars with log r c g > 3.79, we take the lower parts of the isochrones (below the turn-off points). 
A comparisons of values of surface gravity thus derived with those adopted in the high resolution spectro- 
scopic analyses is given in Fig. |l(c)| 

Comparing to values provided in SDSS DR8, our estimates of the surface gravity show much less scatter 
with those adopted in the high-resolution spectroscopic analyses. 

2.3.3 Microturbulent Velocity 

Given the limited resolution and S/N of the SDSS spectra, we derive the microturbulent velocity by its 
dependence on effective temperature and surface gravity. Guided by values adopted in the high-resolution 
spectroscopic analyses, we determine the microturbulence values using Equation (9) in Edvardsson ( 1993) 
for log T c s (K) > 3.79. For log T e g (K) < 3.79, a constant microturbulent velocity v tur b = 1.5 km s _1 is 
adopted. 
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Fig. 1 (a) Comparison of Ti and Ti effective temperatures deduced from the (g — z)o color and those adopted in 
high-resolution analyses for the 13 metal-poor sample stars. The effective temperatures given in SDSS DR8 are plotted 
in (b). (c) Comparison of surface gravities log g obtained from the isochrones and those adopted in high-resolution 
analyses, (d) Same as (b) but for surface gravities log g. 



2.3.4 Model Atmospheres 

We utilize the local thermodynamical equilibrium (LTE) and line-blanketed Kurucz model atmospherej^j 
By interpolation in the discrete grids of model atmospheres, we obtain model atmospheres corresponding to 
the effective temperatures and surface gravities derived above for our sample stars. We employ the program 
ABONTEST8 for abundance analysis. ABONTEST8 assumes LTE, and includes effects of natural broadening, 



thermal broadening, van der Waals damping and microturbulent broadening (Zhang & Zhao 2005 1 



2.4 Derivation of Metallicity 

Using the photometric temperatures and the EW of the line as input, and [Fe/H] from DR8 as initial metal- 
licity, we use the ABONTEST8 program to calculate the initial estimate of the a-element abundance. The 
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Fig. 2 The 12Gyr Yonsei-Yale isochrones! Demarque et al. 2004) for the 13 metal-poor sample stars, adopting the 
metallicities given in SDSS DR8 as initial metallicities. The isochrones (solid lines) are in order of increasing metallicity 
from left ([Fe/H] = —3.3) to right ([Fe/H] = —2.1). The dashed lines show the value of surface gravity corresponding 
to the effective temperature estimated for each star. For cooler stars with log T c g < 3.79, we take the upper parts of 
the isochrones (above the turn-off points), while for hotter stars with log 77 e ff > 3.79, we take the lower parts of the 
isochrones (below the turn-off points). 

iron abundance is then calculated assuming [a/Fe] = +0.4. This new iron abundance is then used to recal- 
culate the isochrone, the surface gravity and model atmosphere. The new set of surface gravity and model 
atmosphere, along with the new iron abundance, are then used as the second set of input parameters to 
ABONTEST8. The process is iterated, until the derived iron abundance differs from the last input value by 
less than 0.1 dex. 

There is no significant difference between the metallicities estimated using the two sets of photometric ef- 
fective temperatures T-y and For the Ca II K line, T 2 provides results that are in slightly better agreement 
with the high-resolution spectroscopic values. Accordingly, we have adopted the metallicity values derived 
fromT 2 . 

2.5 Decision on Final Metallicity Results of the 13 Stars 

Our results show that only the metallicities derived from the Ca II K line agree well with those determined in 
high-resolution spectroscopic analyses. The metallicities estimated from the Mg lb lines show a large scat- 
ter relative to the high-resolution spectroscopic results (see Fig. |3(b)[ ). And the metallicities deduced from 
the Ca II triplet lines are much too high (see Fig. |3(c)| >. The Ca II K line is the strongest metal absorption 
line over the optical wavelength range for most stars, while the Mg I b and the Ca II triplet lines have much 
weaker line strengths. For the metal-poor stars of interest here, given the low resolution and limited S/N of 
the SDSS spectra, these metal lines are difficult to measure due to their weakness (the EWs can be as low as 
~ 100 mA), and small errors in fitting the continuum can lead to large errors in the line strength estimates. 
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The effects are particularly profound for weak features. Apart from measurement uncertainties, NLTE ef- 
fects may also account for parts of the abundance estimate differences, especially for the Ca II triplet lines 
(Mashonkina et al. 2007 ; Starken biirg et al. 2010} . To account for this, we have used the methods presented 
in Gehren et al. (120061 and Mashonkina et al. d20071 > to perform NLTE calculations for the Mg I b lines and 
Ca II triplet lines, respectively. 

Figs. |3(a)[ |3(b)| and |3(c)| compare the metallicities derived respectively from the CallK, Mglb and 
Call triplet lines, without and with the NLTE effect corrections. The abundances determined from the 
Call triplet lines remain higher than the high-resolution values, even after the NLTE corrections. The dis- 
crepancies are probably due to measurement uncertainties. The big circles in the plot indicate the two stars 
that are very enhanced in [Ca/Fe] or [Mg/Fe]. SDSS J234723 has a [Ca/Fe] ratio of +1.11 (ILai et al. 20091 1 
and SDSS1707+58 has [Mg/Fe]= +1.13 ( lAoki et al. 2008K Such high [a/Fe] values can also contribute 
to the discrepancies in the estimated [Fe/H]. Table [2] summarizes the averages and standard deviations of 
the metallicities of our sample derived from each metal line, along with those from SDSS DR7 and DR8. 
Table |2] shows that the iron abundances deduced from the Ca II K line have a higher consistency with the 
high-resolution spectroscopic results than those derived from the other two sets of metal lines. Therefore, 
we have adopted the metallicities based on the measured EWs of the Ca II K line as the final iron abundances 
of our sample. The last three columns of Table [TJ list the effective temperatures (T2), surface gravities and 
iron abundances obtained in this work. 



Table 2: Averages and standard deviations of the differences between the [Fe/H] estimated from the 
different metal lines, those given in DR7 and DR8 and those from the high-resolution spectroscopic 
analysis of the 13 stars. 
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0.4 


0.4 


0.6 


0.5 


0.3 


0.2 


Standard Deviation 


0.3 
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3 SEARCHING FOR EMP STAR CANDIDATES 

As shown above, our method that makes use of the Ca II K line performs well in returning a reliable and 
accurate estimate of metallicity, that is, the method is effective in measuring the metallicities of EMP stars 
using low-resolution spectra. We therefore apply this method to a large sample of metal-poor stars selected 
from the SDSS archive. 

3.1 Selection of the Sample and Measurement of EW Values 

We select 166 EMP stars (including those with repeated observations) from the SDSS archive which have 
[Fe/H] < —3.5, as listed in DR8, and spectra with the average S/N per pixel of more than 20 over the 
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wavelength range 4000 — 8000 A. Hot stars, with effective temperatures higher than 7500 K, are excluded. 
Using the method described above, we derive their metallicities based on measured line EWs. The results 
derived from the weak Mg I b and Ca II triplet lines are not used for the reasons described above. 
For EMP stars with [Fe/H] < —3.5, the line features are very weak. Even for the relatively strong CallK 
line, it is difficult to identify the feature and fit the continuum using low-resolution spectra. Also, it should 
be noted that the assumption of [a/Fe] = +0.4 may lead to overestimated [Fe/H] for stars with strongly 
enhanced a element. Nevertheless, we have managed to measure the EWs of the Ca II K line and determine 
the metallicities for 45 stars out of the 166 EMP candidates. 

3.2 Results and Discussions 

From the results above, we have identified six metal-poor star candidates with [Fe/H] ~ —4.0 using our 
method. It is noteworthy that one of them, SDSS J102915+172927, was recently confirmed to have a metal- 
licity as low as [Fe/H] = —4.73 (based on a ID model atmosphere analysis) using the Very Large Telescope 
(VLT) high-resolution spectroscopic observations (Caffa u et al. 201 il l. The discovery of this low-mass and 
most metal-poor star, without the enrichment of C, N and O, sheds light on the theory of formation of the 
first generation of low-mass stars in the Milky Way. Compared to the VLT result, our iron abundance is 
overestimated by ~ 0.6 dex. The discrepancy is mainly caused by contamination of the stellar Ca II K line 
by interstellar absorption, which cannot be resolved in the low-resolution SDSS spectra. The other EMP star 
candidates with [Fe/H] < —3.8 in our sample remain to be confirmed. Table[3]lists the stellar parameters of 
the six EMP star candidates that we have found. 



Table 3: Parameters of the six EMP Star Candidates. 
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+121152.0 


+010313.1 


-001542.7 


+230803.0 
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Dec. 
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T eB (DR8) 
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5150 


4462 


6129 


6294 
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log g (DR8) 
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[Fe/H] (DR8) 


-3.8 


-3.5 


-3.7 


-3.7 


-3.6 


-3.5 


r e ff (This work) 


5775 
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4435 
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log g (This work) 


3.4 
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4 SUMMARY 



Using the EW of the Ca II K line, effective temperature deduced from the (g — z)o color and surface gravity 
obtained from the isochrone, we can derive a reliable estimate of metallicity for EMP stars with the SDSS 
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spectra. 

We find the metallicities derived from the Ca II K line have good agreement with the high-resolution spec- 
troscopic results. From a sample of 166 potential EMP stars selected from SDSS DR8, we identify six 
metal-poor star candidates with [Fe/H] ~ —4.0, and one of them, SDSS J102915+172927, has recently been 
confirmed to be ultra metal-poor ([Fe/H] < —4.0) by the VLT observations (Ca ffau et al. 201 il l, showing 
our method is effective even at the lowest known metallicities. 

Further high-resolution spectroscopic analyses for the other stars with estimated metallicities [Fe/H] 
~ —4.0 in our sample is of considerable interest. With a rich source of low-resolution spectra provided 
by the ongoing SDSS and the forthcoming LAMOST campaigns, our method has a great potential to dis- 
cover more EMP stars in the near future. 
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Fig. 3 (a), (b), (c) Comparison of our metallicities of the 13 stars derived from the EW of the Ca II K line, Mg I b lines 
and Ca II triplet lines with the values determined by high-resolution spectra, (d) shows the comparison between SDSS 
DR8 and high-resolution results as a reference, (b) and (c) also show our results after NLTE corrections. The big circles 
indicate stars with large over-abundances of Ca or Mg over Fe. 



